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Abstract 


We present an updated catalog of 46,753 radial velocity (RV) standard stars selected from the APOGEE DR17. 
These stars cover the Northern and Southern Hemispheres almost evenly, with 6296 being red giants and 3896 
being main sequence stars. These RV standard stars are stable on a baseline longer than 200 days (with 54% longer 
than one year and 10% longer than five years) with a median stability better than 215 m s^ !. The average number 
of observations of those stars is 5 and each observation is required to have signal-to-noise ratio greater than 
50 and RV measurement error smaller than 500 ms '. Based on the new APOGEE RV standard star catalog, we 
have checked the RV zero-points (RVZPs) for current large-scale stellar spectroscopic surveys including 
RAVE, LAMOST, GALAH and Gaia. By careful analysis, we estimate their mean RVZP to be 4-0.149 km gt 
+4.574 kms! (for LRS), —0.031 kms! and +0.014 km s ^, respectively, for the four surveys. In the RAVE, 
LAMOST (for MRS), GALAH and Gaia surveys, RVZP exhibits a systematic trend with stellar parameters 
(mainly [Fe/H], Terr, log g, Gap — Grp and Grys). The corrections to those small but clear RVZPs are of vital 
importance for these massive spectroscopic surveys in various studies that require extremely high RV accuracies. 
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1. Introduction 


The velocity component of a star in the line of sight direction 
can be defined by the Doppler shift of the spectrum captured by 
the telescope. It can be converted into the framework of the 
solar system's center of mass, called the “barycentric” or 
“heliocentric” radial velocity (RV). The “barycentric” or 
"heliocentric" RV represents the rate of change of the distance 
between the Sun and the star (for a detailed definition of RV, 
see Lindegren & Dravins 2003). The measurement of RV is 
essential to the construction of complete stellar 6D information 
(3D position and 3D velocity). Its accuracy is required to be 
better than several km g or even a few m sl, for various 
Galactic studies such as understanding the structure and 
assembly history of the Milky Way (Binney & Tremaine 
1987; Bovy et al. 2012b; Bovy 2015; Belokurov et al. 2018; 
Gaia Collaboration et al. 2018b; Helmi et al. 2018; Helmi 
2020), estimating the mass of the Milky Way (Xue et al. 2008; 
Bovy et al. 2012a; Huang et al. 2016; Eilers et al. 2019; Zhou 
et al. 2022), defining orbital parameters and characteristics of 
binary systems (Mayor & Queloz 1995; Gao et al. 2017; El- 
Badry et al. 2018; Tian et al. 2018; Liu et al. 2019; Li et al. 
2022b), identification of exoplanets (Xie et al. 2016; Trifonov 
et al. 2020) and systematically searching for hypervelocity stars 
(Brown et al. 2005, 2014; Huang et al. 2017, 2021; Koposov 
et al. 2020; Li et al. 2021, 2022a, 2022c; Marchetti et al. 2022). 


In the past decades, RVs have been measured for over tens 
of millions of stars from a series of large-scale spectroscopic 
surveys, including, ground-based surveys, such as GALactic 
Archaeology with HERMES (GALAH, De Silva et al. 2015; 
Buder et al. 2018), Sloan Digital Sky Survey/Apache Point 
Observatory Galactic Evolution Experiment (SDSS / APOGEE, 
Blanton et al. 2017; Majewski et al. 2017; Abdurro’uf et al. 
2022), RAdial Velocity Experiment (RAVE, Steinmetz et al. 
2006; Kunder et al. 2017), SDSS/Sloan Extension for Galactic 
Understanding and Exploration (SDSS/SEGUE, Yanny et al. 
2009; Eisenstein et al. 2011; Dawson et al. 2013; Rockosi et al. 
2022), Large Sky Area Multi-Object Fiber Spectroscopic 
Telescope (LAMOST, Cui et al. 2012; Deng et al. 2012; Zhao 
et al. 2012; Luo et al. 2015) and space-based surveys, i.e., Gaia 
Radial Velocity Spectrometer (Gaia RVS, Katz et al. 2004; 
Cropper et al. 2018; Gaia Collaboration et al. 2018a, 2018b). In 
the near future, more stellar RVs will be obtained thanks to the 
ongoing/planning massive spectroscopic surveys, such as the 
SDSS-V (Kollmeier et al. 2017), 4-metre Multi-Object 
Spectrograph Telescope (4MOST, de Jong et al. 2019) 
and Dark Energy Spectroscopic Instrument (DESI, DESI 
Collaboration et al. 2016; Dey et al. 2019). 

The measurement of RV can be influenced by various 
factors, including the type of instrument, spectral resolution, 
accuracy of wavelength calibration, methodology used to 
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derive RV, and even observation conditions and environments. 
These factors can lead to significant variations in RV 
measurements. To correct for these effects, it is necessary to 
construct a set of RV standard stars which are stable enough in 
a long observation baseline. 

At present, over tens of thousands of RV standard stars have 
been defined by various efforts, including about 5000 bright 
RV standard stars with extreme stability of 15 ms ^! over an 
average baseline of six years constructed by a long monitoring 
project (Udry et al. 1999a, 1999b; Crifo et al. 2007, 2010; 
Chubak et al. 2012; Soubiran et al. 2013, 2018) and over 
18,000 standard stars with a median stability of 240 m s^! over 
a one-year baseline with a large color and magnitude range 
constructed from the APOGEE Data Release (DR)14 (Huang 
et al. 2018, hereafter H18). However, the number spatial 
density of current RV standard stars is so low that it is hard to 
calibrate the RV zero-points (RVZPs) of the RV measurements 
from future massive spectroscopic surveys. 

This paper is an update to H18. Thanks to the long-term 
repeated observations and more southern stars observed during 
SDSS-IV, the number of RV standard stars has been trebled 
with a much larger sky coverage from APOGEE DRI7, 
compared to the previous version of H18. The paper is 
structured as follows. In Section 2, we correct the possible 
RVZPs of the RV measurements of APOGEE DRI17. In 
Section 3, we describe the details of selections of RV standard 
stars from APOGEE DR17. In Section 4, we use the selected 
APOGEE RV standard stars to calibrate the RVZPs of the 
RAVE, GALAH, LAMOST and Gaia surveys. Finally, we 
conclude in Section 5. 


2. Corrections to APOGEE RV Measurements 


As found in H18, the measured RVs of APOGEE surveys 
exhibit a systematic trend as a function of Tere. To correct this 
trend, 1611 reference RV standard stars, collected from the 
various literatures (Udry et al. 1999a, 1999b; Nidever et al. 
2002; Chubak et al. 2012; Soubiran et al. 2013), are adopted to 
calibrate the RVZP of APOGEE DR14 (Abolfathi et al. 2018). 
In this paper, we apply these reference RV standard stars to 
check the RVZP of APOGEE DR17 (Abdurro'uf et al. 2022). 
The details of the compilation of these reference RV standard 
stars are described in H18. Generally, the RVs of these 1611 
reference stars are required to have stability better than 
100 m s^! over a baseline of at least one year. 

By cross-matching the 1611 reference RV standard stars to 
APOGEE DR17, 118 common stars are found to check the 
RVZPs of APOGEE DRI7. The RV differences between 
APOGEE and reference RV standard stars show a significant 
systematic trend along stellar effective temperature (see 
Figure 1). To describe this trend, a simple linear fit is adopted: 


ARV = —1.2146 + 0.2885 x (Tg/10? K) kms-! (1) 
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Figure 1. The differences between the APOGEE RV (RVApogzg) and the RV 
standard star (RVpgr) collected from the literature as a function of effective 
temperature (Tett). Stars with different colors are reference RV standard stars 
selected from different catalogs, as marked in the bottom-left corner. The 
dashed red line is the best linear-fit to the data points. 


The coefficients found here are similar to those reported in 
H18, implying the robustness of the instruments and RV 
measurements. 

We also cross-matched a set of 4813 RV standard stars 
constructed by Soubiran et al. (2018) with APOGEE DR17 and 
205 common stars were left. The systematic trend found by 
these common stars is generally consistent with that displayed 
in Figure 1 for stars with T«r > 4000 K. For cold stars with 
Ty; < 4000 K, few stars are found. 


3. APOGEE Radial Velocity Standard Stars 
3.1. APOGEE Survey 


APOGEE (Zasowski et al. 2013, 2017; Majewski et al. 
2016; Majewski et al. 2017) is a large-scale high-resolution 
(R ~ 22,500) spectroscopic survey in the near-infrared (H- 
band 1.51-1.70 um), provided by the 2.5 m Sloan Foundation 
Telescope (Gunn et al. 2006) and the 1 m New Mexico State 
University (NMSU) Telescope (Holtzman et al. 2010) at 
Apache Point Observatory (APO) in the Northern Hemi- 
sphere, and the 2.5 m Irénée du Pont Telescope (Bowen & 
Vaughan 1973) at Las Campanas Observatory (LCO) in the 
Southern Hemisphere. The APOGEE survey is an important 
part of the SDSS-III (Eisenstein et al. 2011) and SDSS-IV 
(Blanton et al. 2017) programs. The APOGEE survey is called 
“APOGEE” or “APOGEE-1” in SDSS-IIL, and it is called 
“APOGEE-2” in SDSS-IV. APOGEE-1 started its data 
collection in 2011 and ended it in 2014. The SDSS DRIO 
publicly released the APOGEE-1 three-year data set, which 
was subsequently followed by two additional releases in 2015 
and 2016. This accomplishment successfully fulfilled the 
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Figure 2. Hammer projection in R.A. and decl. of the APOGEE RV standard stars. The black line signifies the Galactic plane, on which we mark the positions of the 
Galactic center (GC) and the Galactic anti-center (GAC). The positions of M31, Large Magellanic Cloud (LMC) and Sagittarius Dwarf Spheroidal Galaxy (Sgr dSph) 


are also marked. 


stated objective of observing over 100,000 stars with a 
limiting magnitude of H — 12.2 mag and spectral signal-to- 
noise ratio (S/N) greater than 100. APOGEE- is a 
constituent program of the SDSS-IV initiative, which 
commenced in 2014 and finished in 2021. In addition to 
collecting the data in the Northern Hemisphere, APOGEE-2 
also adopted data from the 2.5 m Irénée du Pont Telescope 
mounted on LCO to expand the observation sky coverage to 
the Southern Hemisphere. The recently released SDSS DR17 
(Abdurro’uf et al. 2022) incorporates the latest version of the 
APOGEE survey, including more than 657,000 stars, and this 
version is also the final version of all APOGEE-1 and 
APOGEE-2 data. The measurement of RV has an uncertainty 
of l00 ms ' and a zero-point offset of 500 ms ^! (Nidever 
et al. 2015). Typical uncertainties for Teg, log g and [Fe/H] 
are better than 150K, 0.2dex and 0.1 dex, respectively 
(Mészáros et al. 2013; García Pérez et al. 2016). 


3.2. Selecting RV Standard Stars from the APO- 
GEE DRI7 


Following H18, to select RV standard stars from APOGEE 
DR17, we defined the weighted mean RV (RV), internal error 
of RV (erv), RV weighted standard deviation (ay) and 
uncertainty of RV (ogy) for each star separately: 


1. RV = Mee, where w; is the weight assigned by the 
i=1Wi 


individual RV measurement error e;, that is, 1/ e and n 
is the total number of observations; 


em 
2. lay = seit: 

ERV Shw , 
3. oky = —— 37 w(RV — RV); 


Oo = 
RV ew? — SW 


4. egy = max(onv/ VN, Ipv/VN). 


We utilize the symbol AT to denote the time baseline and MJD 
to represent the mean Modified Julian Date of the n 
observations. In the following step, we select RV standard 
stars based on the following criteria: AT > 200 days, n 2:3, 
S/Niow > 50 and egy < 200 m s |, where S /Niow represents the 
lowest S/N for the multiple spectroscopic visits to each star. 
Through the above cuts, a total of 46,753 APOGEE RV 
standard stars were selected. The spatial distribution is depicted 
in Figure 2, with full sky coverage. The distributions of time 
baseline AT, number of observations n and ogy for these RV 
standard stars are plotted in Figure 3. Their AT are all greater 
than 200 days (with 5446 longer than one year and 1046 longer 
than five years). The average number of observations for these 


stars is 5. The median opy of all standard stars is 71.75 ms, 


corresponding to a median stability (30gy) of 215.25 ms !, 
better than 240 m s^! of the H18 sample. We show the color- 
(absolute) magnitude distributions of these stars, see Figure 4 
for H against J — K,, and Figure 5 for absolute Mg against 
Ggp — Grp. Due to the selection effect of the APOGEE survey, 
81% of RV standard stars are redder than 0.5 in terms of 
J— K,. Figure 5 contains 30,268 RV standard stars with 
measurable distances, all of which have been corrected for 
interstellar extinction using the 2D dust map from Schlegel 
et al (1998). We employ the empirical relation 
Mg = 3.53 x (Ggp — Grp) — 0.06 to distinguish giants from 
main sequence dwarf stars (see dashed line) Among them, 
62% are red giants and 38% are main sequence dwarf stars. We 
list 46,753 APOGEE RV standard stars including their name, 
H, J — Kj, Tar, RV (after RVZP correction by Equation (1)), 
lerv, Orv, n, CO gy, AT and mean MJD information in Table 1. 
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Figure 3. The time baseline distribution (left panel), observation number distribution (middle panel) and RV weighted standard deviation (ogy) distribution (right 


panel) of 46,753 APOGEE RV standard stars. 
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Figure 4. The color J — K,-magnitude H diagram of APOGEE RV standard 
stars. The top and right insets are histogram distributions along the J — K, and 
H axes, respectively. 


4. Calibrations of Radial Velocity Scales for Large- 
scale Stellar Spectroscopic Surveys 


Next, we use these 46,753 APOGEE RV standard stars to 
check the RVZPs of the RAVE, LAMOST, GALAH and Gaia 
surveys. The calibration results are shown in Figures 6, 7 and 
Table 2. 

(i) The RAVE survey: The RAVE survey has collected 
520,781 medium-resolution spectra (MRS, R ~ 7500) centered 
on the Cal triplet (8410-8795 A) range. The survey has 
released 457,588 individual stars randomly selected from the 
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Figure 5. The color Ggp — Grp-magnitude Mg distribution of APOGEE RV 
standard stars. The black dashed line is the empirical line to separate giants and 
main sequence dwarf stars. 


Southern Hemisphere stars with 9<J< 12 using the multi- 
object spectrograph 6dF on the Australian Astronomical 
Observatory’s 1.2m UK Schmidt Telescope. Estimations of 
RV, atmospheric parameters (Tor, log g and [Fe/H]) and a 
element abundances were described in Kunder et al. (2017). 
We cross-matched RAVE DRS5 with our 46,753 APOGEE 
RV standard stars, resulting in a total of 1284 common stars 
with S/N>10. The comparisons show a mean ARV 
(APOGEE RVs minus RAVE) of +0.149 kms !, with a 
standard deviation of 1.358 km s '. We display the systematic 
trends of ARV with Ter, log g, [Fe/H] and S/N in Figure 6. 
There is no obvious systematic trend of ARV with Ty, log g or 
S/N, however, ARV exhibits a weak linear trend with [Fe/H]. 
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Table 1 
The Final Sample of APOGEE Radial Velocity Standard Stars 

Name H J-K, Tete RV Terv ORV n ORV AT Mean MJD 

(K) (km s^ P) (km s^ )) (km s ^P) ( kms )) (days) (—50000) 
J00:00:00.20-19:24:49.9 10.74 0.40 5501 18.698 0.0223 0.0926 3 0.0535 1091 7904 
J00:00:00.324-57:37:10.3 10.64 0.42 6162 —21.057 0.0261 0.1424 6 0.0581 716 6328 
J00:00:00.68+57:10:23.4 10.13 0.65 5031 —12.914 0.0165 0.1284 3 0.0742 716 6389 
J00:00:05.07+56:56:35.3 8.72 0.77 4981 4.348 0.0114 0.1125 3 0.0503 421 8651 
J00:00:05.35+15:04:34.4 11.17 0.59 4914 18.229 0.0149 0.0275 3 0.0159 2493 7416 
J00:00:12.11-19:03:38.3 10.53 0.36 5740 23.782 0.0250 0.0727 3 0.0420 1091 7904 
J00:00:12.17-19:49:30.6 10.77 0.39 5645 —5.845 0.0192 0.0427 3 0.0247 1091 7904 
J00:00:12.434-55:24:39.1 11.48 0.84 4659 —114.291 0.0156 0.0367 3 0.0212 2954 6813 
J00:00:13.62-19:13:04.2 11.06 0.37 5555 —30.755 0.0223 0.0585 3 0.0338 1091 7904 
J00:00:20.05--57:03:46.8 9.79 1.06 4355 —52.964 0.0080 0.0789 8 0.0279 2972 7296 
J17:15:48.00-30:18:23.9 8.73 1.60 3392 —55.954 0.0150 0.0793 3 0.0458 327 8493 
J17:15:48.77+59:10:10.1 11.05 0.60 4832 —18.619 0.0207 0.0211 4 0.0106 844 6958 
J17:15:49.98+42:38:09.7 10.85 0.61 4887 11.785 0.0170 0.0514 3 0.0297 324 6143 
J17:15:50.95-42:57:40.5 10.18 1.16 4596 —51.058 0.0176 0.0176 3 0.0102 355 8582 
J17:15:51.314-23:55:24.4 10.26 0.59 4949 19.250 0.0162 0.0544 3 0.0314 314 6521 
J17:15:51.77+30:59:06.2 10.80 0.65 4708 —61.426 0.0098 0.1100 8 0.0389 1358 7726 
J17:15:52.42+65:42:59.0 10.19 0.20 6162 —21.094 0.0266 0.1295 4 0.0648 391 8469 
J17:15:52.48+56:43:37.6 9.15 0.74 4565 4.726 0.008 1 0.1234 11 0.0372 1183 8446 
J17:15:52.74+29:07:36.8 11.25 0.84 4455 —81.656 0.0221 0.1205 4 0.0603 455 7856 
J17:15:52.99+30:52:49.0 8.90 0.23 6151 —6.513 0.0152 0.1330 10 0.0420 1359 7740 


This trend can be described by ARV = 0.1058 + 0.4175 x 
[Fe/H] (see Table 2). 

(ii) The LAMOST survey: LAMOST is a 4 m quasi-meridian 
reflecting Schmidt telescope (Cui et al. 2012). The telescope is 
equipped with 4000 fibers distributed in a field of view with a 
diameter of 5°. Within one exposure, LAMOST can obtain 
4000 optical low-resolution spectra (LRS; R~ 2000; with 
wavelength coverage between 3700 and 9000 A) or MRS 
(R ~ 7500; with two wavelength windows of 4950-5350 A and 
6300-6800 A, respectively). 

Over ten million LRS have been released in the recent DR9 
of LAMOST (http://www.lamost.org/dr9/v1.1/). A total of 
7,060,436 stars in the AFGK Stellar Parameters catalog of 
LAMOST DR9 LRS have measurements of RV and stellar 
atmospheric parameters, which are derived by the official 
stellar parameter pipeline: LAMOST Stellar Parameter Pipeline 
(LASP; Luo et al. 2015). To check the RVZPs of the LAMOST 
LRS RVs, we cross-matched the LAMOST DR9 AFGK Stellar 
Parameters catalog with the RV standard stars; 15,600 common 
stars are found with S/N > 10 (average S/N = 96, see Table 
2). The comparisons show a mean ARV (APOGEE RVs minus 
LAMOST) of +4.574 km s land a scatter of 3.844 km s™!. No 
obvious systematic trends with Tor, log g, [Fe/H] or S/N are 
detected for LAMOST LRS RVs (see middle panels of 
Figure 6). 

The MRS parameter catalog released in LAMOST DR9 
contains measurements of stellar atmospheric parameters and 
RV for over 1.6 million stars from 8 million MRS. To check 


the RVZPs of the LAMOST MRS RVs, we cross-matched the 
LAMOST DR9 MRS parameter catalog with the RV standard 
stars, resulting in 6431 common stars with S/N > 10 (see Table 
2). By comparing their RVs (measurements from LASP) with 
the standard stars, multiple peaks are found in the RV 
difference distribution. There are two dominant main peaks, 
one occurring before 2018 October 19 (MJD — 58,410) and the 
other after this date. Prior to 2018 October 19, the mean ARV 
(APOGEE RVs minus LAMOST) was 6.813kms ' with a 
standard deviation of 1.202 km s 1, while after that date, the 
mean ARV was 0.727kms ! with a standard deviation of 
1.183 km s^ '. The main reason for such a significant transition 
in mean ARV arises from the use of different wavelength 
calibration lamps. Prior to 2018 October 19, the Sc lamp was 
employed to calibrate the wavelengths of the LAMOST test 
observation spectra, whereas the Th-Ar lamp has been used 
since then (Wang et al. 2019; Zhang et al. 2021). LAMOST 
MRS provide zero-point corrected RV measurements, with the 
aforementioned offsets largely corrected. If considering the 
formal survey started from 2018 October 19, the offset- 
corrected RVs from LAMOST MRS agree very well with those 
of the APOGEE RV standard stars, with a nil zero-point and a 
scatter of 1.05 kms_' (see middle panel of Figure 6). However, 
the mean RV differences still show a systematic trend with Tae 
(see middle panels of Figure 6). This trend can be described by 
a fourth-order polynomial in Tsp (see Table 2). 

(iii) The GALAH survey: The GALAH survey is a large-scale 
stellar spectroscopic survey. The aim is to collect high- 
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Figure 6. RV differences of common stars (gray dots) between APOGEE RV standard stars and different surveys (top panel: RAVE; middle two panels: LAMOST 
LRS and MRS (zero-point corrected RV measurements after 2018 October 19); bottom panel: GALAH) as a function of T. (first column), log g (second column), 
[Fe/H] (third column) and S/N (fourth column). The black dots and their error bars in each subpanel signify the mean and standard deviation of ARV, respectively, in 
bins of individual parameters (Terr, log g, [Fe/H] and S/N). Blue dashed lines represent the mean difference as shown in Table 2. Red lines trace polynomial fits of the 
RV difference as functions of stellar parameters (the coefficients are given in Table 2). 


resolution (R — 28,000) spectra of approximately one million 
stars in the optical band (four discrete optical wavelength 
ranges: 4713-4903, 5648-5873, 6478-6737 and 7585-7887 À) 
using the HERMES spectrograph installed on the 3.9 m Anglo- 
Australian Telescope (AAT) at the Siding Spring Observatory 
(De Silva et al. 2015). In the third data release (DR3, Buder 
et al. 2021), GALAH provided a total of 678,423 spectra of 
588,571 unique stars, including measurements of RVs, stellar 
atmospheric parameters and individual element abundances. 
We cross-matched the APOGEE RV standard stars with 
GALAH DR3 to examine the RVZP of the GALAH RVs. A 
total of 1839 common stars with S/N > 10 were found, with 
average S/N of 40 (Table 2). The mean value and standard 
deviation of the ARV (APOGEE RVs minus GALAH) are 
—0.031 kms ! and 0.299kms !, respectively. Figure 6 
(bottom panel) shows the systematic trends of ARV with Teff, 
log g, [Fe/H] and S/N. It can be seen from the plot that ARV 
has no trend with [Fe/H] or S/N. However, ARV exhibited a 
curved trend with both T. and log g. Through our validation, 
we find that the systematic trend of GALAH RVs is dominated 


by logg. The trend can be described by a sixth-order 
polynomial about log g, the coefficients of which are presented 
in Table 2. 

(iv) The Gaia survey: The European Space Agency (ESA) 
satellite Gaia (Gaia Collaboration et al. 2016) recently released 
the DR3 (Gaia Collaboration et al. 2022), which provides 
astrometric and photometric data for more than 1.8 billion 
sources. Compared with Gaia DR2, Gaia DR3 provides more 
than 33 million stars with measurements of RV (Katz et al. 
2022) and more than 470 million stars with measurements of 
atmospheric parameters (Fouesneau et al. 2022). The median 
value of RV measurement accuracy is 1.3 kms™! at Ggvs = 
12 mag and 6.4kms™! at Grvs = 14 mag. The RVZP of the 
Gaia DR2 has a systematic trend with Grys, which shows 
ARV —0 kms ' at Grys = 11 mag and ARV = 0.40 kms! 
at Grys = 14 mag (Huang et al. 2018; Katz et al. 2022). 

To check the RVZP of the Gaia DR3, we cross-matched our 
APOGEE RV standard stars with Gaia DR3 to obtain 43,214 
common stars with 3200 K < Tert < 6400 K and S/N > 5. The 
comparison shows a tiny offset of +0.014 km s™! (APOGEE 
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Figure 7. Top panel: Similar to Figure 6, but this figure presents common stars between APOGEE RV standard stars and Gaia DR3, subplots of which show a 
systematic trend of ARV with Grys (left), Ggp — Grp (middle) and Nops (right). Red lines trace polynomial fits of ARV with Grys or Ggp — Gpp, and the coefficients 
are listed in Table 2. ARV’ in the middle panel signifies the RV differences after the corrections to Grys dependent systematics. Cyan line is a second-order 
polynomial taken from Katz et al. (2022) to describe the systematic trend of ARV with Grys when Grys > 11 mag. Bottom panel: The RV precision of Gaia DR3, as 
a function of Ggvs, is reported through comparison with APOGEE RV standard stars, corresponding to metal-rich stars ([Fe/H] > —0.1, left), middle-level metallicity 
stars (—0.4 < [Fe/H] < —0.1, middle) and relatively metal-poor stars ([Fe/H] < —0.4, right). The different colored curves are calculated for each effective 
temperature range of the APOGEE RV standard stars as labeled in the top-left corner on each panel. 
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Table 2 


Comparisons of RVs yielded by the APOGEE RV Standard Stars and RAVE, LAMOST, GALAH and Gaia Surveys 


ARV 
(km s~') 


+0.149 
+4.574 


—0.006 


—0.031 


+0.014 


s.d. (km s^!) 


1.358 
3.844 


1.053 


0.299 


0.561 


S/N 
56 
96 


61 


40 


N 


1284 
15600 


6431 


1839 


43214 


Calibration Relationship 


ARV = 40.1058 + 0.4175 x [Fe/H], for —1 < [Fe/H] < 0.5 
ARV = +4.57 
ARV = —02233 — 0.1772 x [(Teş — 3900)/10?] 


+ 1.3234 x [(Tag — 3900)/107]? — 1.2789x 

KT — 3900)/10?P + 0.3138 x [CT — 3900)/10?]^, for 3900 < Tete < 6500 
ARV = —0.18344-1.3138 x log g — 

1.9784 x log g? + 1.3765 x log g? — 0.4979 x log g4+ 

0.0879 x log g5 — 0.0059 x log g°, for 0.8 < log g < 4.7 
ARV = 0.0163 — 0.0335 x (Grvs ox: 8.2) + 0.1327 x (Gnvs arx 82» 

— 0.0930 x (Grvs — 8.2) 
+0.0232 x (Grys — 8.2)* — 0.0020 x (Grys — 8.2), for 8.2 < Grys < 14 
ARV" = —0.7751 + 1.0987 x (Gap = Grp) — 0.4549 x (Gap = Grp)? 


+ 0.0450 x (Gap — Grp)? 
+0.0064 x (Gap — Grp)’, for 0.5 < Gpp — Grp < 3.3 


è Here ARV’ denotes the RV differences after the corrections to Grys dependent systematics. 
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RVs minus Gaia), with a small scatter of 0.561 kms !. The 
systematic trend of ARV with color Ggp— Grp, magnitude 
Ggvys and number of transits (Nop;) is depicted in Figure 7. 
Significant systematic trends for ARV with color and 
magnitude are detected. For Ggp — Grp, systematic deviations 
are clearly detected at Ggp — Grp < 1 mag and Ggp — Grp > 
2 mag. For Grys, a systematic trend is significantly found at 
Grys > 11 mag. We first adopt a fourth-order polynomial to 
correct the trend along with Grys. After corrections to Grys 
dependent systematics, the trend along with Ggp— Grp is 
further corrected by a fourth-order polynomial fit. The resulting 
coefficients are presented in Table 2. It is worth noting that 
Katz et al. (2022) also identified the Gg ys dependent trend and 
provided a second-order polynomial to describe it for 
Grys > 11 mag (as shown in Figure 7). The second-order 
polynomial can partially capture the systematic trend as we 
discovered, but it cannot correct the systematic trend at the faint 
end (Ggys > 13 mag). 

Based on the common stars of our APOGEE RV standard 
star and Gaia DR3, we study the precision of Gaia RV 
measurements. As displayed in the bottom panels of Figure 7, 
the precision of Gaia RV measurement generally decreases 
with Tert and Grys and slightly increases with [Fe/H]. The RV 
precision in the bright range (Ggys « 10 mag) is several 
hundred ms ', and is a few kms! at the faint end (Grys > 
12 mag). This result is consistent with the prediction of Katz 
et al. (2022). 


5. Summary 


We have constructed a catalog of 46,753 RV standard stars 
from the 657,000 near-infrared (H-band; 1.51—1.70 um) high- 
resolution (R ~ 22,500) spectra provided by APOGEE DR17. 
They are almost evenly distributed in the Northern and 
Southern Hemispheres, with 62% red giants, and 38% main 
sequence dwarf stars. They were observed with a time baseline 
of at least 200 days (with 54% longer than one year and 10% 
longer than five years) and were observed more than three 
times. The median RV stability was 215.25 ms '. Using the 
catalog of RV standard stars, we calibrated the RVZPs of four 
large-scale stellar spectroscopic surveys: RAVE, LAMOST, 
GALAH and Gaia. By careful comparisons, we found the mean 
RVZPs are +0.149 kms !, +4.574 kms ! (for LRS), 
—0.031 kms! and +0.014 kms !, for RAVE, LAMOST, 
GALAH and Gaia, respectively. In addition to an overall 
constant offset, RVZPs of part of these surveys show moderate 
dependences on stellar parameters (e.g., Ter, logg, [Fe/H], 
color or magnitude). We further provide corrections by simple 
polynomial fits with coefficients listed in Table 2. Our studies 
demonstrate that the small but clear RVZPs in these large-scale 
spectroscopic surveys can be well detected and properly 
corrected by our RV standard stars, which are believed to be 
useful for their further applications in various investigations. 
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The complete APOGEE RV standard star catalog in Table 1 is 
publicly available at https:/ /nadc.china-vo.org/res/r101244/. 
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